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Abstract —Model Predictive Control (MPC) can efficiently 
control constrained systems in real-time applications. MPC 
feedback law for a linear system with linear inequality con¬ 
straints can be explicitly computed off-line, which results in an 
off-line partition of the state space into non-overlapped convex 
regions, with affine control laws associated to each region of 
the partition. An actual Implementation of this explicit MPC in 
low cost micro-controllers requires the data to be “quantized”, 
i.e. represented with a small number of memory bits. An 
aggressive quantization decreases the number of bits and the 
controller manufacturing costs, and may increase the speed 
of the controller, but reduces accuracy of the control input 
computation. We derive upper bounds for the absolute error 
in the control depending on the number of quantization bits 
and system parameters. The bounds can be used to determine 
how many quantization bits are needed in order to guarantee 
a specific level of accuracy in the control input. 

I. INTRODUCTION 

Model predictive control (MPC) [1] is an efficient method 
for control design of multivariable constrained systems in 
chemical and process control, automotive, aerospace, and 
factory automation [2]-[4]. MPC solves a constrained op¬ 
timal control problem in real time (on-line). 

Explicit MPC (EMPC) [5], [6] may reduce on-line com¬ 
putational costs and code complexity by pre-computing the 
MPC feedback law as a state feedback, thus making it viable 
for fast applications with limited computational capabili¬ 
ties [7]-[10]. In particular, for linear systems subject to linear 
constraints and cost function based on 1-norm, c»-norm, or 
squared 2-norm, the EMPC results in a polyhedral piecewise 
affine (PWA) feedback law. Thus, during the on-line execu¬ 
tion, the EMPC controller first identifies which polyhedral 
region contains the current state, and then computes the 
control action by evaluating the corresponding affine control 
law. The identification of the polyhedral region is referred to 
as the point location problem [11], which can be solved by 
sequential search and binary search tree see, e.g., [12], [13]. 
Due to the exponential increase of the number of regions with 
respect to the number of constraints in the MPC problems, 
techniques for reducing complexity of the EMPC feedback 
law while maintaining its most important properties have 
been proposed, see, e.g., [14]-[16] and references therein. 

In practice, the data of EMPC have to be typically stored 
in a micro-controller hardware memory, so that every stored 
number is represented by a small fixed number of bits for 
every number in the data. In other words, the data cannot 


be stored exactly and hence a precision loss occurs. We call 
this reduction of precision “quantization” and the reduced 
precision data “quantized” data. The method for quantiza¬ 
tion can be as simple as rounding. Aggressive quantization 
has the advantage of decreasing memory requirements and 
increasing the speed of the control input evaluation, at the 
price of introducing inaccuracy in the computation of the 
control input. If the quantization precision is too small, the 
controller can fail to accurately determine the region for the 
current state of the controlled system, and thus, the control. 
Eor example, by quantizing the state measurement/estimate 
data the quantized state may jump to a different region. 

The effect of quantization has been investigated for im¬ 
plicit MPC for instance in [17], [18]. We investigate the 
resulting accuracy in the control input computation in EMPC 
as a consequence of different quantization precisions, so that 
we can determine how many bits need to be used to guarantee 
a desired level of accuracy in the control input. A brief 
overview of our approach and results is in Section 

In Section we provide a mathematical accuracy anal¬ 
ysis, depending on mutual positions of the exact and the 
quantized system states. When the quantization does not 
affect the system state region, so that the same feedback law 
applies to both the exact and the quantized system states, an 
error bound is easy to establish. A difficult case for analysis, 
leading to a much larger possible controller inaccuracy, is 
where the quantization makes the system state to jump over 
a region facet to a different region. In this case, bounding 
the accuracy of the control requires taking into account not 
only quantization precision for the system state, but also 
quantization effects of the region facets and of the feedback 
laws in different regions. We derive two kinds of upper 
bounds on the accuracy of the control input computation in 
Section]^ Bounds without knowledge of the quantized data, 
describing the worst case scenario, for that reason are called 
“a priori.” After a quantized implementation of the controller 
is determined, the a priori bounds are improved, using the 
already known off-line quantized data, in addition to the 
original data. The resulting tighter “a posteriori bounds” 
depend on the quantization precision of the current state. 
We show how our bounds can be improved by exploiting a 
rescaling technique that makes the system state space evenly 
sized in all spacial directions. 


We validate the bounds numerically in Section IV 


^Mitsubishi Electric Research Laboratories; 201 Broadway Cambridge, 
MA 02139, USA 

^knyazev@merl.com 
®dicairano@merl.com 
Accepted to IEEE CDC 2015 conference 


II. Brief overview 

Throughout this paper, R denotes the set of real numbers, 
II • 111, II • II 2 , and II • ||oo denote 1-norm, 2-norm, and oo-norm, 
respectively. A' denotes the transpose of A. 



The control u{x) = {ui(a;),... ,Un^{x)} is a continuous 
PWA function determined by the EMPC control law 


symbol “hat” to denote the data after the quantization, and 
the symbol “A” for the quantization error. 


— ^iX Gi Vx G Pit t — 1, . . . , (1) 


where a; € K", the gains Fi G and offsets Gi G K™, 

and n-r denotes the number of regions 


P, = fxG U’^lffiX < Ku H, G 


with iJ, = 


-1 ! 


\K,GW^'^}t 




and K, = [Kl..., ]'. 


The controller determines on-line, i.e. in real time, which 
region contains the given state x. This is typically the most 
time consuming operation, for large n and rir, requiring 
computing numerous matrix-vector products HiX. If the state 
x is in the region Pi, then the corresponding control law ([T]) 
is used to compute the control input u{x). All the controller 
“true” data, determining the regions, and control law 0 are 
computed off-line, typically in the double precision computer 
arithmetic, and then are quantized and stored in a memory 
of the controller as quantized data. 

The reduction of the precision of the data used by the 
controller decreases the number of bits stored by the con¬ 
troller and increases the speed of the control, but reduces the 
accuracy of the on-line execution of the controller. The target 
accuracy of the controller can be based on the accuracy of 
the sensor for sensing, or of the estimator for estimating the 
curi'ent state of the system. The state of the controlled system 
determined with quantized data deviates from the state of the 
controlled system determined with the true data within the 
limits depending on the control law, the data representing 
the system, and the precisions of the quantization. Thus, 
for any reduction of the precision of the quantized data, 
it is possible to determine off-line bounds for a maximal 
deviation of the state of the controlled system caused by 
that reduction. Using these bounds, different reductions in 
the precision, compared to the true precision, can be tested, 
and the maximal reduction of the quantization satisfying the 
accuracy requirement of the control can be selected. 

It may also be advantageous to use different quantization 
precisions for various data. For example, the EMPC control 
law Ui{x) = FiX + Gi is evaluated only once on-line, for 
the determined index i, so the gains Fi and offsets Gi can 
be quantized in high precision and stored in a slow-access 
memory without noticeably affecting the controller speed. 
Moreover, different facets can benefit from using different 
numbers of bits in their quantized format to represent the 
same level of accuracy in the control, while speeding up the 
computationally challenging point location on-line search. 

However, in some bounds below, for simplicity of pre¬ 
sentation we suppose that all data are quantized using the 
same scalar quantization function f{z), assuming that z = 
f{z) = z + Az, where |Az| < e and 0 < e < 1, aiming 
at a fixed point quantization, rather than a floating point 
rounding. From the true state x, we obtain the quantized 
state X = X + Ax, and, similarly, Ui{x) = FiX + Gi, 
where HiX < Ki for i = 1,..., We commonly use the 


III. Accuracy analysis 

In this section, we focus on accuracy analysis of the 
control input computation. The control input depends on the 
location of the state vector, thus, to analyze the accuracy of 
the control input computation we have to examine how the 
state changes and which region it falls into before and after 
the quantization. When it falls into the region with the same 
index, i.e. x is in the region Pi and x is in the region Pi, 
the absolute accuracy of the control input is measured by the 
maximum absolute changes between Ui{x) and Ui (x), i.e., 
||ui (x) — iti(x)||oo, and can be easily bounded above by the 
precisions of the quantization. 

The control computation error is generally larger if the 
original state and the quantized state belong to different 
regions with large and different gains. We analyze the case 
where x is in the region Pi and x is in the region Pj, with 
i ^ j, which means the state x is in one of the regions, 
but after quantization the state x gets into another region, so 
the accuracy of the control input is ||ui (x) — itj(x)||oo, as 
illustrated in Fig. [1] where n = to = 1. Fig. shows that 
the change in the slope (gain) when the state jumps over the 
hyperplane (a point in Fig. [T] , since n = 1) can increase the 
control computation error. 



Fig. 1. Location of the state within the regions. Top: the true data, where 
the true state x is in the true region Pi giving the true control u. Bottom: the 
quantized data, where the quantized state x jumps into the quantized region 
Pj having a different gain(slope) of the quantized control u, resulting in a 
large error u — u. 


However, the jump can only happen if a distance from the 
state vector to the facet separating the two regions is small, 
and the PWA control function is continuous, which allows 
us to bound the error even in the case of the jump. 

The sharp upper bound of ||u(/(x)) -u(x)||^ is 

max max \\Fi f{x) + Gi — FjX — Gj\\oo- 

xGPj,f{x)GPi 

We derive an analytic bound for a parr i and j under 
simplifying assumptions, e.g., that Pi and Pj share a facet. 









We start with analyzing the accuracy of the hyperplane 
representation, where the data are quantized. 

Lemma 3.1: Let hx < k he a half-space. Let y = hx — k 
and y = {h + Ah){x + Ax) — {k + Ak) with ||A/i||oo < L 
||Aa;|joo < e, and |Afc| < e for some e > 0. We have 

\y-y\<S := e{\\h\\i + ||a;||i + ne + 1). (2) 

Proof: By direct calculation, we obtain 

\y — y\ = \{h + Ah){x + Ax) — {k + Ak) — {hx — k)\ 

= \hAx + Ahx -f AhAx — Ak\ 

< \hAx\ + |Aft,x| + \AhAx\ + |AA:| 

< ||/i||i||Aa;||oo+ 

||A/i||oo||a:||i -f ||A/i||oo||Aa:||i -f |Afc| 

< e||/i|ji + e||a;||i + ne^ + e. 


Corollary 3.2: Let the hyperplane hx — k border two 
neighboring regions Pi and Pj , such that hx < k for x G Pj 
and hx > k for X G Pi. Let the state x G Pj, but after the 
quantization X G Pi with i f j. Then —S <y<0<y<S. 

Proof: By Lemma \y — y\ A After quantization, 
X falls into the region Pi which means x is out of region Pj. 
So, 6 + y > 0. Therefore, we have —S < y <0. Similarly, 
we can obtain that 0 <y <S. m 

Corollary 3.2 states that, as a result of the quantization, 
the state can jumps to the other side of a hyperplane only if 
the distance from the state to the hyperplane is less than 6, 
as illustrated in Fig. |^for n = 2. 



Fig. 2. The state x in Pi after quantization turns into x in Pj. 

We are now ready to prove our main result. 

Theorem 3.3: Let the hyperplane hx = k border two 
neighboring regions Pi and Pj, such that hx < k if x G Pj 
and hx > k if X G Pi. Let the state x be in the region Pj, 
the orthogonal projection of x on the hyperplane hx = k he 


in a facet of Pi, and x = f{x) G Pi. Then 

\\u{x) - u(a;)||oo < pjll- Fj)h'\\oo + 

\\FiAx + AFiX + AFiAx + AGj||oo, 
where 5 = e(||/i||i + ||a;||i + ne + 1) and further 
\\FiAx + AFiX + AFiAx + AGi||oo 

— II AFi||oo||a;||oo + ||AGi||oo + ||.Gi||ooe, (3) 

or, alternatively, 

\\FiAx + AFiX + AF^Ax F AGj||oo 

^ £(||Ti||oo + ?^||2;||oo + ne + 1). (4) 

Proof: By the triangle inequality, we have 

||M(i) - u(a;)||oo = ||ni(^) - nj(a:)||oo 

= \\FiX + Gi — FjX — Gjjloo 
< \\FiX + Gi — FjX — Gjj|oo+ 

IjF^Aa; + AFiX F AF^Ax F AGi||oo- 

We first bound above IjF^x + Gi — FjX — GjHoo, where x 
satisfies —S < hx—k < 0 by Corollary |3.2| Let Xp denote the 
orthogonal projection of x on the hyperplane hx = k, then 
Xp = xFth', where h' is the transpose of h and f is a scalar. 
We have k = hxp = hx F thh' = hx F t\\h\\2. Therefore, 
|f| = \hx — A:|/|j/i|| 2 . Since u{x) is a linear continuous 
affine function, we have FiXp F Gi = FjXp F Gj, which is 
equivalent to Fi{x F th') F Gi = Fj{x F th') F Gj. Hence, 

\\F,x + G, - F,x - G.lloo = ||^(^’^ - T^O^'lloo 

= |f||l(i^.-F,)/l'||oo 

^ * IIl||2 *II(^*-^j>'II“ 

The second term in the triangle inequality is 

||FjAx + AFiX F AFiAx F AGj||oo 

= ||(Ai^,x + AG,) + (i^*Ax + Ai^,Ax)||oo 

< ||AFiX F AGilloo + II.Fill 00 II Ax ||00 

< IIAFi||oo||a;||oo + ||AGi||oo + ||Fi||ooe. 

Alternatively, we bound the second term as follows 

\\FiAx F AFiX F AF^Ax F AGi||oo 

< ||FjAx||oo + llAFix||oo + ||AFjAx||oo + ||AGi||oo 

< ||Fi||oo||Ax||oo + ||AJ^i||oo||a::||oo 

+ II ALi||oo|| Ax||oo + ||AGi||oo 

< e{||Fi||oo + n\\x\\oc FntF 1). 

Combining the inequalities above completes the proof ■ 
Bound (|^ uses the data obtained before and after quantiza¬ 
tion, therefore we call it a posteriori. In contrast, bound Q, 
called a priori, is more pessimistic, being based on the true 
data only, without the quantized data explicitly appearing. 










In bounds ([^ and Q and in the expression for <5, the 
norms of the state vector can be further bounded, e.g., since 
X € Pj, by the largest norm of any vector in the region Pj. If 
it is not known which region the state belongs to, we can use 
as a universal upper bound the largest norm of any vector in 
the state space, assuming that the latter is bounded. 

Theorem |3.3| implies that the control error is significantly 
reduced if the gains Fi and offsets Gi are quantized with 
high precision. Theorem 3.3 also suggests choosing differ¬ 
ent precisions for quantizing different facets maintaining a 
uniform level of control accuracy over the state space, e.g., 
if ||(Fi — Fj)h!\\oo is small, the hyperplane hx = k can be 
quantized with low precision, without significantly affecting 
the control computation accuracy. 


A. Rescaling 

The accuracy of u{x) depends on how small e and 6 are. 
To make e and S small, one approach is to rescale the regions, 
h, and k, such that 

hD~^Dx k 

max(||hD“i||i, |fc|) ~ max(||/iZ)“i||i, |fc|) ’ 

where is a diagonal matrix with ||I3a;||oo < 1- Let Dx, 
/iD“^/max (||/iD“^||i, |fc|), and fc/max |fc|) 

be our new x, h, and k, respectively. We have ||a;||oo < 1, 
\hx\ < 1, and |fc| < 1. Our new <5 is less than e{n + 2 + ne). 
The control law has the form Ui{x) = FiD~^x + Gi, where 
X & Pi, i = 1,... ,nr. Consequently, 

||?i(i) - u(a;)||oo < Pin“ Fj)D-^h!\\^+ 

^\\FiD + nei\\x\\oo + neei + El, 


where ||Aa;||oo < e = 2~^ for a vector x. The quantization 
errors AHi, AKi, AFi, and AGi are known, given the 
number b of bits for the fractional part. Every component 
in AiJi, AKi, AFi, AGi is bounded by 2~^. 

We generate a consistent uniformly distributed random 
state X, and for every such a state compute the a priori and a 
posteriori bounds of absolute errors of u(x) as described 
in Theorem |3.3| as well as the actual absolute error of 
u{x) obtained by a hxed point model of the controller. 
We eliminate the states with very small a posteriori bounds 
and actual errors (due to the fact that in some regions the 
gain is nearly zero), order the remaining states by sorting the 
a priori errors, and plot in Fig. for two values of a and b. 
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where e is taken, after rescaling, for the regions and ei is 
taken, after rescaling, for the input control. 

Rescaling of the state is evidently equivalent to changing 
the units of the components of the state vector. The bounds 
show that it may be beneficial to choose the units in such a 
way that the state space is balanced in size in all components. 


IV. Test Results 

In this section, we present tests for a double integrator 


x{k + 1) = Ax{k) + Bu{k) 

y{k) = Gx{k) Du{k), 


where 


A = 


1 1 
0 1 


,B = 


,C = 


The state x satisfies the constraints 


■ 1 0 

,D = 

0 ■ 


0 1 

0 


■ -15 ■ 

< X < 

■ 15 

* 

-15 


15 



and the input u satisfies the constraints —1 < m < 1. 

In our tests, we use a fixed point number format which has 
a specific number of bits reserved for the integer part and a 
specific number of bits reserved for the fractional part. We 
use the MATLAB function fi with a-bit total word length, 
1-bit for sign and 6-bit fraction length, such that, e.g.. 


X = fi(x, 1, a, 6) = X + Ax, 


Bounds and absolute errors of u(x) with a=16 and b=9 
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Fig. 3. Evaluating a-bit word total length and b-bit fraction length. Top: 
a priori bounds, a posteriori bounds and real computation absolute errors 
with a = 12 and 6 = 5. Bottom: both bounds of absolute errors of u{x) 
and absolute eiTors in real computation with a = 16 and 6 = 9. 


In Fig. top (bottom) panel for a = 12 (a = 16) and 
6 = 5 (6 = 9), the maximal difference between our a priori 
bounds and actual errors is about 0.9 (0.12), the maximal 
























difference between our a posteriori bounds and actual errors 
is about 0.26 (0.03), and the maximal actual error is about 
0.1 (0.02), i.e., approximately 300% (1000%) compared to 
the quantization error 2“^ « 0.03 (2“® « 0.002). 

We observe in Fig. bottom panel that the bounds and 
the actual errors form clusters. These clusters correspond 
to different hyperplanes, demonstrating that it may be ad¬ 
vantageous to quantize hyperplanes using a plurality of 
precisions aiming at a uniform behavior of the bounds and 
the actual errors of the control over the state space. Next, 
we specifically test the states near the common hyperplanes 
between some pairs of neighbor regions to classify these 
hyperplanes in terms of their sensitivity to the quantization. 

The results are summarized in Table For example, the 
bounds and the real errors are less than 10 “^, if the state 
is in region Pi or Pg and after quantization the state jumps 
into region Pg or Pi. Checking the controller data, we find 
that the gains in the regions Pi and Pg are nearly zero, thus, 
naturally, the control does not change if the state is in Pi and 
Pg, which is also well captured by our a posteriori bound, 
thus the hyperplane separating the regions Pi and Pg can be 
quantized with very low precision. We notice several other 
trivial pairs of the regions in Table |I] identified by the small 
bounds and errors. Let us check closely one nontrivial pair. 


Neighboring 

regions 

Maximal 

a posteriori bound 

Maximal 
actual error 

1 and 9 

less than 10~^ 

less than 10”^ 

2 and 8 

less than 10“^ 

less than 10“^ 

3 and 5 

0.19 

0.15 

3 and 12 

0.25 

0.12 

4 and 7 

less than 10”'‘ 

less than 10“^ 

5 and 6 

less than 10~^ 

less than 10”^ 

6 and 9 

less than 10“^ 

less than 10“^ 

7 and 8 

less than 10“'‘ 

less than 10“^ 

7 and 13 

0.21 

0.12 

10 and 13 

0.26 

0.13 


TABLE I 

A POSTERIORI BOUNDS AND REAL ERRORS FOR THE STATES NEAR A 
SINGLE COMMON HYPERPLANE WITH a = 12 AND 6 = 5 

e.g., P 3 and P 5 . The maximal a posteriori bound is 0.19 and 
the largest actual error is 0.15, if the state is in region P 3 or 
P 5 and after quantization the state jumps into region P 3 or 
P 5 . Fig. 1^ plots the a posteriori bounds and the actual errors 
for this case. We observe that our a posteriori bounds always 
bound above the actual errors and are reasonably sharp. 

Another nontrivial example is the P 3 and Pi 2 pair. The 
maximal a posteriori bound is 0.25 and the largest actual 
error is 0.12. Fig. plots the a posteriori bounds and the 
actual errors for the P 3 and P 12 pair. Fig. plots the similar 
data for the P 7 and P 13 pair, and Fig. for the Pig and P 13 
pair. We observe in these figures that our a posteriori bounds 
overestimate the actual errors about two times. 

We finally note that our bounds in all tests depend on the 
norm of the state, not on the state itself Bounding above the 
norm of the state would make the bounds state-independent. 


Bounds and absolute errors of u(x) with a=12 and b=5 



Fig. 4. A posteriori bounds and real errors with a = 12 and 6 = 5 for 
the states near the common hyperplane between regions 3 and 5. 


Bounds and absolute errors of u(x) with a=12 and b=5 



Fig. 5. A posteriori bounds and real errors with a = 12 and 6 = 5 for 
the states near the common hyperplane between regions 3 and 12. 

Conclusions 

EMPC data are quantized and stored in the memory of 
the controller. A state-of-the-art method for determining 
the precision of the operations in the controller uses a 
subjective decision based on an ad hoc educated guess of an 
engineer designing the controller for the given system. Such a 
subjective decision does not usually guarantee a specific level 
of accuracy of the control. The validity of the quantized data, 
representing the system in the controller is typically checked 
numerically versus the true data on randomly selected state 
vectors. The state space in true precision includes so many 
vectors that an exhaustive validation is impractical even off¬ 
line. Inaccurate on-line computation of the control can result 
in suboptimal control, system malfunctions, and failures. 


























Bounds and absolute errors of u(x) with a=12 and b=5 



Fig. 6. A posteriori bounds and real errors with a = 12 and b = 5 for 
the states near the common hyperplane between regions 7 and 13. 


Bounds and absolute errors of u(x) with a=12 and b=5 



Fig. 7. A posteriori bounds and real errors with a = 12 and 6 = 5 for 
the states near the common hyperplane between regions 10 and 13. 


The actual EMPC control errors can be much larger 
compared to the quantization error, if the state is near a 
hyperplane separating neighboring regions with large and 
different gains, requiring a special analysis. We propose 
restricting random state vectors to the neighborhoods of the 
region facets, dramatically decreasing the off-line computa¬ 
tional time for numerically checking the control accuracy 
versus the true data. We analyze the EMPC control accuracy 
deriving upper bounds under an assumption on the quantiza¬ 
tion error, typical for a fixed point arithmetic, commonly used 
in controllers to improve on-line performance. An influence 
of a rescaling of the state space on the accuracy of the 
control computation is examined. Using our bounds, one 
can determine the required precision of the quantization and 


estimate the accuracy in the control input, designing the 
controller. It is shown that various EMPC data have different 
sensitivity with respect to the precision of the quantization, 
e.g., the gains and offsets are preferTed to be stored with 
high precision. We discover that it may be advantageous 
to use different precisions quantizing various hyperplanes, 
where the data larger affecting the accuracy of the controller 
are stored with a greater precision, compared to the data 
less affecting the accuracy of the controller. Numerical tests 
for a simple double integrator system support our theory 
and conclusions. Our future work concerns probabilistic 
approaches, taking into account that the state vector may 
follow trajectories rather than filling the whole state space. 
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